Field-Induced Ising Criticality and Incommensurability in Anisotropic Spin-1 Chains 
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Quantum phase transitions induced by an external magnetic field in the Haldane-gapped spin-1 
chains are studied in a fermionic field-theoretic description of the model. In the case with broken 
axial symmetry, two transitions occurs consecutively as increasing the field: h c = (mim^) 1 / 2 (with 
mi and mi the transverse masses) is the critical Ising point and h* (> h c ) is the commensurate- 
incommensurate transition point. Although the latter is thermodynamically indiscernible, its prop- 
erty is revealed in the spin correlation functions (both uniform and staggered) which are calculated 
upon expressing the correlations of coupled Ising pair fields in terms of block- Toeplitz determinants. 

PACS numbers: 75.10.Pq, 75.40.-s, 75.30.Gw, 71.10.Pm 



It is now well established that a one-dimensional (ID) 
quantum antiferromagnet with spin-1 possesses a (Hal- 
dane) gap in the excitation spectrum which separates the 
first excited triplet states from the singlet ground state 
0. The elementary excitations are coherent antiferro- 
magnetic massive magnons appearing around q = tt and 
carrying quantum number S = l. If an external magnetic 
field is introduced into the Haldane system, the trans- 
verse modes couple to the field and phase transition may 
occur when the field overcomes the gap 0, D IE EI- 
Experimentally, there have indeed been observed phase 
transitions in some quasi-lD compounds by measurin; 
thermodynamic quantities such as magnetization 00 
and specific heat |3, llfij . as well as neutron scattering 
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For an isotropic chain or the field applied having ax- 
ial symmetry, a commensurate-incommensurate (C-IC) 
transition [ljj takes place when the field surpasses the 
transverse mass, and the system enters the massless Lut- 
tinger liquid phase. Physically, this is understood in ci- 
ther bosonic 0,0 or fermionic language 0. In the former 
description, the transition is associated with the conden- 
sation of the bosons; while in the latter, the fermions 
begin to fill the conduction band when the field (as chem- 
ical potential) exceeds the gap. Aside from magnetiza- 
tion, the external field induces a planar quasi-long-range 
order in the staggered spin correlations. The correlation 
functions in the Luttinger liquid phase have also been 
calculated |f|. However, in most real materials the spin 
rotational symmetry is hardly preserved in presence of 
single-ion anisotropies, and thus the triplet masses are at 
least slightly different. Considering the Luttinger liquid 
requires a U(l) symmetry and cannot survive a broken 
axial symmetry, in a general asymmetric case, the phase 
transition reduces to be of Ising- type . In this Letter, 
we report the results for the Ising scenario, in particular, 
the spin correlation functions in all phases. 

We take a field-theoretic model due to Tsvelik 
for the anisotropic Haldane chain: ~Hs(x) — 

E a=1 , 2 ,3[-f ~ &d*t.Z) ~ imaggf], where v 

is the characteristic velocity of magnons, and the triplet 



£u = (£i,<^,£^) are chiral {v = R,L) Majorana (real) 
fermion fields. The masses, m a (> 0), are consid- 
ered as phenomenological parameters accounting for the 
anisotropy. Although the neglected marginal current- 
current interactions are essential to give rise to field de- 
pendent exponents in the Luttinger liquid phase, their 
role is exhausted in the Ising transition case by renormal- 
izing the masses and velocity. Since the system is now 
massive everywhere except right at the critical point, the 
effect of interactions is at most shift the transition point 
without changing the universal property significantly. In 
this representation, the spin density, as a sum of uniform 
and staggered parts, is expressed as S = I + (—l) n jN 
(7 being a constant). The uniform part I = Ir + II 
with the chiral currents J„ = — ^^X^ satisfying the 
SU(2)2 Kac-Moody algebra and having the scaling di- 
mensions (0, 1)/(1, 0) for v = R/L. The staggered part is 
associated with the primary field with scaling dimension 
(j|,j|): N = (0-1/12^3,^10-2^3, M1M203), where a a and 
fi a are, respectively, order and disorder Ising variables. 

The magnetic field coupling to the uniform spin den- 
sity contributes a Zeeman energy as Hm(i) = hl z (x) = 
—ihi^x^ii + ClCl), which apparently only mixes first 
two (transverse) species of Majoranas, and the third one 
drops out of the system as a free massive Majorana, or 
equivalently, off-critical Ising model. Therefore, we need 
only concentrate on the Hamiltonian 



a=l,2 



(1) 



Denoting m = |(mi + 7712) and A = |(mi — 1112) and 
combining two Majoranas into a Dirac fermion, model 
JTJ actually describes a charge density wave-Cooper 
pairing-chemical potential (m-A-h) system in ID. When 
bosonized, it becomes a sine-Gordon model with extra 
terms of dual field and field derivative. 

While the uniform spin correlations are easy to com- 
pute, it is not the case for the staggered part, since they 
are composed of Ising order and disorder fields, which are 
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related to the Majorana fermions in a highly non-trivial 
way ( Jordan- Wigner transformation) . To exploit the full 
advantage of the solvableness of the model, we regularize 
on a ID lattice, on which the relations between the 
Ising variables and the fermions can be unambiguously 
established. By using the relations between the Majo- 
rana fields and the lattice Dirac fermion: 



3_(e i7r /4.4 



0, 



+ ^( e ~ W4 4 + e" /4 c„). 



(2) 



we have the following lattice version of two coupled quan- 
tum Ising chains (QICs) (or m-A-h model): 

H = W + H W _ ih y ( a l K _ 6 t an) , (3) 



where H, 



QIC 



_ ./ 



En [(4 _ c «)(4+l + c «+l) - A c( c 



c ")(4 + c n)] with parameters J = ^ and A c = 1 + ^ 
(c= (1, 2) = (a, b)). Without loss of generality, we assume 
m\ > m 2 > (so that m > A > and A > 1) and h > 
unless otherwise stated. The units are so chosen that 
J=a=l. 

Model © is diagonalized through a unitary transfor- 
mation and the spectrum is found to be 



where 8 2 = 



E± = (e 2 ± 2Ti 



-A 2 + h 2 and r 2 : 



,1/2 



-M) 1/2 



(4) 



with 



A| = e 2 + sin 2 fc and ej, = cos k — \. There exist totally six 
(real) correlation functions for the fermions: 

Cf(n) = (cico), C 2 » = <44)> ( c = a,b) 
C 3 (n) = iia^bo) , C 4 (n) = -i(a n b ) . (5) 

The explicit forms of these correlations are obtained in 
a variational approach, and all of them are ID integra- 
tions over momentum. To calculate the correlations of 
the coupled Ising fields, we need further 



Ga,i 



-6 n0 + 2[C^ b (n)-C^' b (n)}, (real) 



F n = 2i[C 3 (n) - C 4 (n)] 



(imaginary) (6) 



As is well-known, the correlation functions of the QIC 
can be expressed as the so-called Toeplitz determinants. 
A Toeplitz matrix has its elements parallel to the diago- 
nal the same value, i.e., it takes the form T= [fr(i — ])]■ 
We find that for the two QICs coupled in the present 
fashion, the correlation of a pair of Ising fields can be 
cast into a block- Toeplitz form with a doubled size. To 
this purpose, we introduce following seven rxr Toeplitz 
matrices: 



[-Fi- j+1 ], F ll (r) = [F i - j ], 

-F-i+j-l] > 



(7) 



where the indices (i, j) run from 1 to r, and the matrix el- 
ements G°' b and F n are given by Eqs. 10 • Then, the Ising 
pair-field correlations like C® b (r) = (cri^ 2 (n)aiyU 2 (n + 0) 
reduce to (2r) x (2r) block- Toeplitz determinants: 



CfJr) 



det 
det 
det 
det 



G%(r) -F a (r) 
F,(r) [Gl(r)] T 

G b a {r) -F a {r) 
[F,(r)f [G«(r)] T 

G%(r) -F„(r) 
F,(r) [Gl(r)] T 

6%{r) -F„(r) 
F,(r) [G b „{r)\ T 



(8) 



With all these provisions, the spin correlations are 
readily obtained. For the uniform part, Sf{r) = 
{I a (n)I a (n + r)): 

S*> y (r) = S rO [C b - a (Q)+C c 1 (0)] 

~2C b 1 ' a {r)C c 1 {r)-2C h 2 a {r)C c 2 {r), 
Sf (r) = 5ro[Ct{0) 

-2[C 3 (r)] 2 +4[C 3 (0)] 2 + 2[C4(r)] 2 . (9) 

Here Cf 2 ( r ) are f°r the third (decoupled, off-critical) 
QIC. The extra terms in Sf(r) is due to the admixture of 
the first two QICs in the presence of the magnetic field. 
Actually, 4[C 3 (0)] 2 is the magnetization squared. For the 
staggered part, S%(r) = (N a (n)N a (n + r)>: 



S x N (r)= C% (r) C; (r), S* N (r) = Cf a (r) C£ (r) , 



S z N {r)=CZ{r)C c a {r) 



(10) 



where C c a „ (r) are Ising correlations for the third decou- 
pled QIC, which bear simple Toeplitz forms (rxr matri- 
ces). 

The fundamental and instructive quantity is the en- 
ergy gap. The gap of the first two coupled QICs (m- 
A-h model) is determined by Ag b p = mm(E^). In ax- 
ially symmetric case (A = 0), the gap decreases lin- 
early: Ag b p (h) — (m — h)6(m — h). In general case, 
the gap function is shown in Fig. ^ (solid line), which 
has the following features: When h is small, it di- 
minishes from m 2 (the smaller mass of the two) as 
Agap(' 1 ) R:;m 2 — 2^s- I n the vicinity of the critical point, 
h c = (m 2 — A 2 ) 1 / 2 = (mim 2 ) 1 / 2 , the gap vanishes linearly 
at both sides as A^ (h) « ^\h — h c \. When the mag- 



netic field is further increased to h* 
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position in momentum space begins to shift from k = to 
KW « [{ 2 -f ) 1/2 + (h ~ K) 1 ' 2 . As a result, 

a gapful incommensurate phase is expected at large h. 
On the other hand, it is shown that when h > h* and 
for small A, the gap Ag^ p ~A. Therefore, in the axially 
symmetric case, the gap closes in this region. 



FIG. 1: The energy gap (A^ p ) of the coupled QICs (m-A-h 
model) as a function of h (the solid line). The gap vanishes 
at h c = 0.0866025 for m = 0.1 and A = 0.05. Above h* c 
(=0.108080), the dashed line is for the wavevector at which 
the gap locates. The (numerically extracted) solid circles rep- 
resent the effective mass, m e «, that governs all correlation 
functions; and the open circles are for the true incommensu- 
rate wavevector, k*, that enters the correlation functions in 
the incommensurate phase. The vertical dash-dotted line at 
h* = 0.098664 is the true C-IC boundary. The wedge (dotted 
line) stuck to the h c point represents the mass deduced from 
the effective Ising theory. 

The identification of the Ising criticality at h c fol- 
lows from the linearly vanishing behavior of the gap and 
the diagonalizcd Hamiltonian. Since E£ always remains 
massive, the modes that associated drop effectively out 
of the low-energy theory, which is resultingly an Ising 
theory keeping only half of the original degrees of free- 
dom. By expanding the spectrum around the transition 

1 /2 

point and for long wavelengths, « J e s (k 2 + m 2 s ) , 
where J e ff= ^ < 1 is a new (smaller) energy scale for the 
effective theory with an effective mass m c g- = ^ | h — h c \ 
(the dotted wedge connected to the Ising point in Fig-QJ. 

The emergent T — phase diagram is shown in Fig. [2 
in which we release the constraint on the relative val- 
ues and signs of the masses, so that it covers the whole 




FIG. 2: Phase diagram of the coupled QICs (m-A-h model). 
The sheets enclosed by the solid lines belong to the Ising crit- 
icality, and those by the dashed lines are the C-IC transition 
sheets. The triangle with dash-dotted sides stands for the 
incommensurate U(l) critical plane (Luttinger liquid). 




FIG. 3: Four Ising pair correlation functions at various mag- 
netic field strengths. The parameters are same as in Fig. Q 

parameter ranges for the m-A-h model {Q or The 
critical Ising sheets, determined by the condition h = 
h c (mi, 7712), arise from the trivial Ising lines (mi = and 
777,2 = at h = 0), and merge the C-IC sheets, pinned 
down by h = h*(mi, 7712)1 at the U(l) symmetric line 
h = mi = 7772, i-e., the usual C-IC transition line in the 
axially symmetric case (A = 0). Above the transition 
only the system with the U(l) symmetry is critical. 

Unfortunately, although A°^ p and k* are enlightening 
and analytically expressible, neither of them is the true 
quantity that controls the asymptotic behaviors of the 
correlation functions. We have to numerically extract the 
corresponding effective mass m c ff (i.e., inverse of the cor- 
relation length) and effective incommensurate wavevec- 
tor k* from the fermion correlation functions, Eqs. J5J, 
by evaluating the integrals. For comparison, we show 
these two quantities in the same Fig. ^ (the solid and 
open circles). As a result, the real C-IC transition point 
h* locates between h c and h* . We also see from the figure 
that the effective mass, which is in general greater than 
the gap of the system, deviates significantly from the gap 
above the Ising transition. The mass grows rapidly un- 
til the system enters the incommensurate phase where 
the mass becomes field independent. In the vicinity of 
the Ising transition, the linearized effective mass (dotted 
wedge) asymptotically coincides with the numerical mass 
values (solid circles). The dissociation between the gap 
and mass is in fact due to the external field which breaks 
the Lorentz invariance of the model. 

By combining Eqs. JHJ, J7J, and 10, and the known 
values of functions Cf(r)-C4(r), we compute the four 
correlation functions of the Ising pairs numerically. The 
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TABLE I: Asymptotics of the uniform and staggered spin correlations when A > 0. 
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results are shown in Fig. EI We find that below the Ising 
transition C^(r), as in the h = case, is asymptot- 
ically a constant, while above the transition C^(r) be- 
comes a constant. The lines decrease faster (in fact twice) 
in C^ h a (r) before the Ising transition and in C^(r) after 
the transition. Furthermore, even in the incommensurate 
phase, functions (r) and C^ a (r) are not modulated by 
the wavevector k*. We also recognize that with respect 
to the effective Ising theory, the combined operators /Ui<72 
and take effectively the roles of the Ising order and 
disorder fields across the transition, since their correla- 
tions have the standard power and exponential forms. 

Finally, by Eqs. © and (|10|l . we obtain the results 
of spin correlations in various phases. In the axially 
symmetric case (A = 0), when the field exceeds the 
C-IC transition point, we find the transverse uniform 
spin correlations: Sj' v (r) cx r~ 3 / 2 e~ m3r sin[fc c r — S(k c )] 
with the phase shift 8(k c ) ~ — and the incommensurate 
wavevector k c r* (^j 1 ) 1 ^^ — m) 1 / 2 , and the longitudi- 
nal uniform spin correlations: ASj (r) — Sj (r) — M oc 
— r~ 2 [l — cos(2fc c r)] with M=^- the magnetization; and 
the transverse staggered correlations: S^ v (r) cx r -1 / 2 , 
and the longitudinal staggered correlations: S^(r) cx 
r 1 e _m3r cos(fc c r). These results are in agreement with 
those from other theories just above the transition 
point. In anisotropic case, the asymptotics for the spin 
correlations are summarized in Table [I] Here, the phase 
shifts in the incommensurate phase, restricted within 
(0, 7r), are in general field and mass dependent; rj{< 1) 
is also a non- universal constant. As argued by Affleck 
0, we verify that one of the transverse staggered spin 
correlations (S V N in our case) behaves simply like the cor- 
relation of an Ising variable across the transition, and the 
fluctuation of the longitudinal uniform correlation obeys 
the power-law A6>f cxr -2 at the criticality. 

To summarize, we have studied the Haldane spin-1 sys- 
tem with generic anisotropies under the effect of external 
magnetic field. We showed explicitly the Ising critical- 
ity in both the spectrum and correlation functions. The 
latter was calculated based on the establishment of the 
block- Toeplitz determinant representation for the corre- 
lations of the Ising pair-fields. A true long-range order is 
formed above the transition. We also found the incom- 
mensurability in high field is a rather robust property, 



which appears in both iso- and anisotropic systems, and 
the mass (not the gap) is field independent in the incom- 
mensurate phase. Finally, we emphasize the fact that the 
C-IC transition (point h*) in the massive phase is ther- 
modynamically unidentifiable, in other words, there is no 
anomaly associated with this transition in either magne- 
tization or specific heat measurement. The discussions 
in this respect, as well as the detailed formalism of the 
present work, will be published elsewhere. 
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